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It is increasingly important to understand, in detail, two-pion correlations measured in p + p and d +A col- 
lisions. In particular, one wishes to understand the femtoscopic correlations, in order to compare to similar 
measurements in heavy ion collisions. However, in the low-multiplicity final states of these systems, global 
conservation laws generate significant A'-body correlations which project onto the two-pion space in non-trivial 
ways and complicate the femtoscopic analysis. We discuss a formalism to calculate and account for these 
correlations in collisions dominated by a single particle species (e.g. pions). We also discuss effects on two- 
particle correlations between non-identical particles, the understanding of which may be important in the study 
of femtoscopic space-time asymmetries. 



I. INTRODUCTION 

The unique and distinguishing feature of collisions between 
heavy ions is their large (relative to the confinement scale) 
size and the possibility to generate bulk systems which may 
be described in thermodynamic terms, allowing to discuss the 
Equation of State of strongly-interacting matter. The pri- 
mary evidence for the creation of bulk matter at the high- 
est energies d S S 01 is the existence of strong collec- 
tive flow @]. The dominant feature of flow is the cor- 
relation between space and momentum which it generates; 
thus, momentum-onlypbservables such as pj spectra and az- 
imuthal anisotropics [ 1, 2, 3, 4] represent only an indirect pro- 
jection of the effect. Femtoscopic measurements access space 
as a function of particle momentum, thus providing the most 
direct probe of the most crucial feature of heavy ion colli- 
sions [c.f. e.g. In particular, flow is manifest by a negative 
correlation between the "HBT radius" and the transverse mass 
(mj) of the particles |7|]. 

Clearly, then, a detailed understanding of femtoscopic mea- 
surements in heavy ion collisions is crucial to proving the ex- 
istence of, or probing the nature of, the bulk system generated 
in the collision. It is in fact possible to quantitatively inter- 
pret both the femtoscopic and momentum-only observations 
at RHIC- in A+A collisions- in consistent, flow-dominated 
models of the system [e.g.|H]. All seems well. 

However, it is important to understand the system size de- 
pendence of the apparent bulk behavior. In this paper we dis- 
cuss the possible complications in the comparison of large and 
small systems. 



Especially since low pj observables directly reflect the bulk 
medium, it is reasonable to ask whether similar comparisons 
between Au+Au and p+p collisions reveal comparable differ- 
ences in the soft sector. 

Common measurements of this type include total particle 
yields, pj spectra and azimuthal correlations iy-f). However, 
in each case, it is not clear if qualitative differences between 
small and large systems are observed. Quantum number con- 
servation laws in small systems can strongly affect particle 
yields [e.g. [I2,[T3|]. However, modulo canonical suppression 
effects, equilibrium-based thermochemical fits to yields from 
p+p collisions produce quite similar results to those from 
A+A collisions 1 1511 . Blast-wave fits (aLLs] to pr spectra from 
p+p collisions indicate sizable radial flow, though smaller than 
that seen in Au+Au collisions lfl7ll . Intrinsic anisotropies from 
p+p collisions are considered as non-flow contribution to az- 
imuthal correlations in A+A. However, it is far from obvious 
that the finite values of V2 from p+p collisions [ 18] do not arise 
from collective flow in the p+p collision itself. 

A Quark-Gluon Plasma is usually considered a form of 
matter. If a Quark-Gluon Plasma were created in p+p col- 
lisions, as suggested by Bjorken lfl9ll . would it display bulk 
properties? A direct comparison of soft-sector observables 
in p+p and A+A collisions is necessary, to address these is- 
sues 0201, l2lll . The imminent hadronic and heavy ion program 
at the LHC brings the relevance of such studies into strong 
relief. 



B. Femtoscopy in p+p collisions 



A. Hadron collisions as a reference to heavy ion collisions 

One of the most exciting signals at RHIC is the modifica- 
tion of the jet structure due to the bulk medium. In particular 
leading particle distributions J^, [Toll and azimuthal correla- 
tions lllll in A+A collisions are strongly suppressed relative 
to those from p+p collisions, at high pj. 
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More light might be brought to bear on this important ques- 
tion through femtoscopic measurements, which probe more 
directly the space-momentum correlations generated in a col- 
lective system. 

Though not as plentiful as in heavy ion collisions, two- 
pion femtoscopic measurements are common in e + + e~ or 
p+p(p) collisions Ir22tl . In these collisions, too, "HBT radii" 
are observed to fall with mj. Speculations of the physics be- 
hind this observation have included Heisenberg uncertainty- 
based arguments, string-breaking phenomena, and tempera- 
ture gradients; an excellent overview may be found in 12311 . 
Preliminary measurements by the STAR Collaboration 112411 
even suggest that p+p collisions show similar collective be- 
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havior as A+A collisions. Distinguishing different physical 
mechanisms, however, requires a detailed understanding of 
the correlations themselves. 

Such an understanding is complicated by the clear obser- 
vation of non-femtoscopic effects in two-pion correlations in 
small systems. For example, in A+A collisions the functional 
form (Gaussian or not) fitted to two-pion correlations[e.g.0] 
incorporates only femtoscopic effects. Such fits for smaller 
systems [e.g. |25[] have required additional ad-hoc terms of 
non-femtoscopic nature. 



C. Non-femtoscopic correlations 

Femtoscopic correlations are those which depend directly 
on the two-particle coordinate-space separation distribu- 
tion [c.f. @]. In general, such correlations are confined to 
low relative velocity. Non-femtoscopic correlations may arise 
from string fragmentation or global conservation laws, for ex- 
ample, and there is no reason to expect that such correlations 
appear only in kinematic regimes (e.g. ranges of relative mo- 
mentum) different than the femtoscopic ones. Thus, separat- 
ing those correlations may be a non-trivial exercise. 

Non-femtoscopic correlations may arise from a variety of 
sources. Jets will clearly induce momentum-space correla- 
tions between their fragmentation products. While these ef- 
fects may not be negligible, the low momentum of the pi- 
ons under consideration {pj ~ 0.4 GeV) puts us squarely 
in the region in which factorization breaks down and the jet 
interpretation becomes significantly murkier. In the kine- 
matic region under consideration, string fragmentation may 
play a role; this is an area for future study, though significant 
model-dependence will be present. Collective bulk flow (e.g. 
anisotropic elliptic flow) will generate A^-body correlations 
which will project onto the two-body space. "Clusters" within 
events- i.e. several independent particle-emitting sources- 
may generate additional structure j26il27ll : indeed, such clus- 
ters may be treated as "large resonances" with a many-body 
decay channel. Each of these sources of non-femtoscopic cor- 
relations may play a greater or lesser role in a collision, de- 
pending upon the physical scenario. In this paper, we do not 
focus on these sources of correlation. 

In this work, we focus on effects which must be at play in 
any physical system - Energy and Momentum Conservation- 
Induced Correlations (EMCICs). These global conserva- 
tion laws provide an A^-body constraint on the event, which 
projects down onto 2-body spaces and should become more 
pronounced at lower event multiplicity (N). 

EMCIC effects on femtoscopic correlation functions have 
been estimated [see Appendix C of |27l] recently in the con- 
text of a numerical model of Bose-Einstein correlations from 
emitting cells, using a rough but fast numerical algorithm to 
conserve energy-momentum. An earlier study by Bertsch 
et al. UM included energy-momentum conservation in an 
analytically-solvable model in the limit of one spatial dimen- 
sion and non-relativistic particles. Both studies were confined 
to correlation functions in one dimension of relative momen- 
tum. In this paper, effects of EMCICs on three-dimensional 



correlation functions are studied in detail, using numerical 
simulations and analytic expressions, both based on the re- 
stricted phase-space integral. 

D. Structure of this paper 

In Section we describe GENBOD , an event generator 
which samples an inclusive momentum distribution subject 
only to constraints of energy and momentum conservation. 
In Section Hill we briefly discuss the harmonic representation 
which provides a complete and natural characterization of the 
shape of the correlation function. An extensive discussion 
of symmetry constraints on the harmonics is given in Ap- 
pendix [B] 

For the next three Sections, we focus on events in which 
only pions are emitted. In Section[IV]we use pion-only events 
from GENBOD to illustrate the effects of varying constraints, 
frames, and kinematic cuts, on EMCICs. A method to cal- 
culate analytically (but using information from measured dis- 
tributions) EMCICs is shown in Section [V] This leads to an 
"experimentalist's formula," given in Section |VI] intended to 
disentangle EMCICs from other (e.g. femtoscopic) correla- 
tions in the data. The formula involves several approxima- 
tions which may break down in reality; these are discussed 
and effects are evaluated quantitatively. 

In Section IVIII we discuss two effects that can compli- 
cate a direct comparison of EMCICs for identical and non- 
identical particle correlations. Several interesting effects are 
observed, which might be important for the increasingly com- 
mon studies of space-time asymmetries. We find that, for non- 
identical particles, the "experimentalist's formula" is only ap- 
proximately applicable. We summarize our discussion in Sec- 
tion [vm] 

IL CALCULATING EVENTS WITH ENERGY AND 
MOMENTUM CONSERVATION 

To clearly understand the role of EMCICs, we would like 
to study events in which there is no other physics involved be- 
sides the conservation laws. Such a tool has been provided 40 
years ago in the form of the GENBOD computer program [see 
l29l for an excellent write-up of the method and physics] in 
the CERN library. Given a requested total multiplicity (AO, 
a list of masses (m,) of emitted particles, and a total amount 
of energy (£ t ot) to distribute among them, GENBOD returns an 
event of random momenta (four-vectors pj), subject only to 
the condition of energy and momentum conservation. More 
importantly, it returns, for each event, a weight proportional 
to the probability that the event will actually occur in nature. 
Thus, it is a much different tool than, say, transport codes like 
RQMD 13011 . in which each event returned may be treated as 
equally probable. 

This weight is based on the phase-space integral Rp/ ll3lll 

r AN ( N \ N 

Rn = J 5 4 (p-£^JnS(p?-m?)A<, (1) 
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FIG. 1: (Color online) A high-probability multiplicity-30 event 
calculated by GENBOD . Lines correspond to particle momenta 

Px,Py,Pz- 

where P = ^E tot , Oj is the total momentum four- vector of the 

event. Rn figures dominantly in Fermi's statistical theory ll32ll . 
in which the probability of having N particles in the final 
state is proportional to Sn ■ Rn', here Sn is the phase-space- 
averaged S-function (or matrix element) associated with the 
process generating the final state. 

In the limit for which momentum distribution is dominated 
by phase-space restrictions alone, S is a constant, and the spec- 
trum of a quant ity a (s ay, an angle or transverse momentum) 
is given by Bill III 

f(a) = ^-R N . (2) 
da 

In the limit that a represents the ensemble of momenta con- 
stituting a given event, Equation |2] returns the event weight. 
See 12911 for a practical iterative prescription to calculate Rn 
and the weights. 

We select (via Monte-Carlo) GENBOD events according to 
their weight and run them through identical software as used 
for experimental analysis. Fortunately, the code is fast, since 
one must calculate large statistics from which to select. This 
is because the phase-space weights vary by large factors. As a 
very extreme case, Figures[T]and|2]show a likely and unlikely 
event, respectively, for multiplicity N = 30. As one would 
expect, the "rounder" event is more likely, though one might 
be surprised by the factor of a hundred million between the 
probabilities. 



III. SPHERICAL HARMONIC DECOMPOSITION OF 
CORRELATION FUNCTIONS 

Measurements of the space-time extent of a particle- 
emitting source at the femtometer scale is commonly done by 
analyzing two-particle correlation functions C (q) as a func- 
tion of the relative momentum q. Experimentally, C (q) is the 



FIG. 2: (Color online) A low-probability multiplicity-30 event calcu- 
lated by GENBOD . Lines correspond to particle momenta p x ,Py,Pz- 

ratio of the q distribution when both particles are measured in 
the same event, to the same distribution when the two particles 
come from different events [c.f. e.g.Q, and references therein, 
for more details]. In this Section, we present correlation func- 
tions produced in exactly the manner in which experimental 
ones are formed. GENBOD -generated events are selected ac- 
cording to the returned event weight. 

Ideally, then, any structure remaining in this ratio reflects 
the correlation between particles in the same event. In the 
present study, these correlations come from energy and mo- 
mentum conservation effects. 

In this paper we will use the commonly used Bertsch-Pratt 
("out-side-long") decomposition of the relative momentum 
q (HHHIl], where q„ is parallel to the transverse total momen- 
tum of the pair, qi is parallel to the beam direction and q s is 
perpendicular to those. 

Usually the 3-D correlation functions are presented in 1- 
D Cartesian projections (or slices) along these axes (e.g. q ) 
with the other q-components (e.g. q s and qi) small. Such 1- 
D slices for GENBOD calculations are presented on Figure [3] 
At asymptotically high relative momentum \q\, femtoscopic 
contributions to the the correlation function (those described 
by the Koonin-Pratt equation [discussed in|^]) must approach 
a constant value, usually normalized to unity, independent of 
the direction of q. Naturally there are no femtoscopic cor- 
relations in these events; correlations induced by the global 
conservation laws are signaled by the non-unity value of the 
correlation function. The correlation function depends on the 
direction of q as well as \q\. 

However, one-dimensional projections represent a set of 
zero measure of the 3-D correlation function and are thus a 
poor tool for exploring its detailed and potentially important 
structure. In principle, one could visualize the full structure of 
the 3-D correlation function via a series of Cartesian projec- 
tions in qt over different ranges in qi^, where i 7^ j 7^ k. This 
would, however, constitute a large number of figures, and rel- 
evant patterns which cut across projections might not stand 
out. 
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FIG. 3: (Color online) ID projections of 3D correlation function cal- 
culated in LCMS frame for multiplicity-9 event calculated by GEN- 
BOD . 



By exploiting symmetries in g-space, the spherical har- 
monic decomposition (SHD) B57I1 becomes a much more ef- 
ficient representation which uses all of the data to show the 
shape of the correlation function. Here, the spherical coordi- 
nates 9, (j), and Q = \q\ relate to the Cartesian ones as 



<?» = 2sin0cos(j), q s = gsinGsincj), qi^Qcosd, 



(3) 



and we define harmonic moments A/ m 's as 
1 



^(cose)c(<2,e ! (|))y/ im (e,(|)). (4) 



Usually, experimentally measured correlation functions are 
not continuous functions of Q, cos 8 and (j), but are constructed 
with bins of finite size. In this case, Equation |4] needs mod- 
ification to account for finite-bin-size effects. For the exper- 
imental practitioner, we discuss one way to deal with this in 
Appendix [C] For the remainder of this manuscript, we will 
assume these binning effects have been dealt with; i.e. we 
assume negligible bin size in cos and (j). 
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FIG. 4: (Color online) Correlation function for the same data pre- 
sented in Fig. [3] is shown at a fixed value of Q = 0.79 GeV/c (ap- 
proximately indicated by the shaded region in Fig. [3} as a function 
of (]) for five bins in cos (9). Curves represent the SHD components 
of various orders; see text for more details. 



Symmetry constrains the number of relevant components. 
For femtoscopic analyses of identical particles at midrapidity 
which integrate over reaction-plane orientation (i.e. almost all 
analyses to date), only real parts of A/ „,'s with even values of 
I and m do not vanish. For the complete list of symmetries of 
A/ m 's, see AppendixlBl Further, it is natural to expect that the 
statistical relevance of high-/ components is diminished. 



As an example, Fig. [4] shows the calculated correlation 
function (the same as shown in Fig. O for one value of Q 
as a function of cos0 and (j). Also shown are curves repre- 
senting SHD with increasingly higher order components. In 
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FIG. 5: (Color online) SHD coefficients for GENBOD -generated 
events consisting of 9 pions having average kinetic energy per par- 
ticle K = 0.9 GeV, as measured in the pair LCMS frame. No kine- 
matic cuts were applied to data. Green squares are A;,,„'s from the 
GENBOD events. Gold inverted triangles are the SHD coefficients of 
Equation [24] for k=2. Black circles, blue stars and red triangles are 
SHD coefficients of the first, second and third terms, respectively, of 
the right side of Equation [25] Open circles are SHD coefficients of 
the right side of Equationl25l 
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FIG. 6: (Color online) SHD coefficients for GENBOD -generated 
events consisting of 18 pions having average kinetic energy per par- 
ticle K = 0.9 GeV, as measured in the pair LCMS frame. No kine- 
matic cuts were applied to data. Green squares are Ai Jn 's from the 
GENBOD events. Gold inverted triangles are the SHD coefficients of 
Equation [24] for k=2. Black circles, blue stars and red triangles are 
SHD coefficients of the first, second and third terms, respectively, of 
the right side of Equation [25] Open circles are SHD coefficients of 
the right side of Equation [25] 



particular, the curves correspond to 

C LM {Q,QA) = V4^( X? E A ljn (Q)Y,: m (Q^) 

\l=0 m=-l 
M \ 

+ E A L , m (Q)YLM) • (5) 

m=—M J 

For example, the curve labeled as "L=2 M=0" contains Ao.o 
(the constant term) and A2.0 components. 

Clearly, for this example, only the first few components are 
required to represent the structure of the correlation function. 
While a few higher-/ terms may be required in some cases, 
the number of relevant A/ m 's is generically expected to be 
small. This is from general considerations of smoothness and, 
for experimental data, statistical issues. Thus, by glancing at 
only a few one- dimensional plots, one views the entire correla- 
tion structure in orthogonal components. The number of plots 
is usually reduced further by symmetry constraints (c.f. Ap- 
pendix [B). 

As an example, the first few A/ „,'s for the same GEN- 
BOD calculations presented in this section, are plotted as a 
function of Q in Fig. [5] The odd-Z and -m moments (not 
shown) vanish as required by symmetry (c.f. Appendix iBli. 



IV. EMCICS FROM GENBOD 

In this Section, we briefly discuss factors which affect the 
A; ,„moments, using Figures IBTTOl For the present, we focus 
only on the green squares, labeled "CF (GenBod)," in those 
Figures. 

Figures [5] [6] and [7] show the A/ ,„'s calculated in LCMS 
frame ^ from GENBOD events that have the same average 
kinetic energy per particle (K = 0.9 GeV) but different mul- 
tiplicity. As expected, the strength of the EMCICs decreases 
with event multiplicity. Similarly, for a given event multi- 
plicity, one expects larger EMCICs when there is less avail- 
able energy. As shown in Figures [6] and [8] for multiplicity- 18 
events, this is indeed the case. 

Since the definition of the "out," "side" and "long" 
directions- and thus the angles 8 and (j)- depend on the frame 
of measurement, one expects the spherical harmonic coeffi- 
cients A/ m to depend on reference frame. This is shown in 
Figures [6] and [9] for correlations measured in LCMS and pair 
CMS frames. 

Less intuitive is the observation that the correlation strength 
depends also on kinematic cuts. Figures l9l and [TOl show the 
A/,„,'s calculated by GENBOD for 18-pion events without and 
with a selection of |t|| < 0.5, respectively. (Note that this cut 
applies to the pions which are used in the analysis, not to the 
set of particles for which energy and momentum is conserved; 
energy and momentum is always conserved for the full event.) 

Finally, we note two important and generic effects. Firstly, 
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FIG. 7: (Color online) SHD coefficients for GENBOD -generated 
events consisting of 6 pions having average kinetic energy per par- 
ticle K = 0.9 GeV, as measured in the pair LCMS frame. No kine- 
matic cuts were applied to data. Green squares are A;,,„'s from the 
GENBOD events. Gold inverted triangles are the SHD coefficients of 
Equation [24] for k=2. Black circles, blue stars and red triangles are 
SHD coefficients of the first, second and third terms, respectively, of 
the right side of Equation [25] Open circles are SHD coefficients of 
the right side of Equationl25l 
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FIG. 8: (Color online) SHD coefficients for GENBOD -generated 
events consisting of 18 pions having average kinetic energy per par- 
ticle K = 0.5 GeV, as measured in the pair LCMS frame. No kine- 
matic cuts were applied to data. Green squares are Ai Jn 's from the 
GENBOD events. Gold inverted triangles are the SHD coefficients of 
Equation [24] for k=2. Black circles, blue stars and red triangles are 
SHD coefficients of the first, second and third terms, respectively, of 
the right side of Equation [25] Open circles are SHD coefficients of 
the right side of Equation [25] 



EMCICs are present at all values of \Q\, reminding us that 
we cannot (responsibly) ignore these effects in a femtoscopic 
analysis. Secondly, in Figures [9] and [TO] we have included 
A; m components up to I = 4. Typically, |A/ + 2, m /A/ „,| ~ 0.1, 
another reminder that characterization of the 3-dimensional 
correlation function requires only a few harmonic compo- 
nents. 



similar to those in Refs IB5L 1361 13711 to obtain correction fac- 
tors which implement EMCICs onto multi-particle distribu- 
tions. In the course of the calculation, we make some simpli- 
fying approximations. The derived expressions are then tested 
for accuracy against the numerical GENBOD simulations. Fi- 
nally, the expressions are used to extract an "experimentalist's 
formula" discussed in the next Section. 



V. ANALYTIC CALCULATION OF EMCICS 

Even if EMCIC effects generated by GENBOD "resemble" 
the experimental data, it is likely unwise to use GENBOD it- 
self to correct the data for several reasons. Firstly, there is 
strong sensitivity to the (not completely measured) number 
and species-mix of all particles emitted in the event, includ- 
ing neutrinos and possible magnetic monopoles (or, less exoti- 
cally, particles escaping detector acceptance). Secondly, there 
is strong sensitivity to the energy "available" in the event; it is 
not obvious that this is y/smf of the collision. Thirdly, EMCIC 
effects depend on the individual momenta p\ and p2 of the 
particles entering the correlation function. This will depend 
on acceptance, efficiency, kinematic cuts and, to a degree, 
the underlying single-particle phase-space. (While correlation 
functions are insensitive to the single-particle phase-space, the 
correlations which they measure may, in fact, depend on this 
phase-space, due to physical effects.) 

Thus, one would like to calculate EMCICs, based on the 
data itself. In this Section, we begin by following arguments 



A. Restricted phase-space corrections 

Danielewicz J35ll . and later Borghini, Dinh and Olli- 
trault l36ll . considered EMCIC-type effects on two-particle 
azimuthal correlations (quantified by V2 and often used as 
a measure of elliptic flow They focused mostly on 

transverse momentum (fir) conservation only, but Borghini 
later l37ll generalized to the case of an arbitrary number D 
of independent (orthogonal) spatial dimensions and recently 
considered momentum conservation effects on three-particle 
analyses of jet-like behavior lf38h . 

As we shall see below, for correlation functions used 
in femtoscopy, the conservation of energy generates ef- 
fects of similar magnitude as those due to conservation of 
(three-)momentum. We deal only with on-shell particles, for 
which energy can/iof be treated as independent of the momen- 
tum (as, say, p x would be largely independent of p y ). Thus, 
unlike the above-mentioned works, we will explicitly begin 
with the more general multivariate central limit theorem. 
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FIG. 9: (Color online) SHD coefficients for GENBOD -generated 
events consisting of 1 8 pions having average kinetic energy per par- 
ticle K = 0.9 GeV, as measured in the pair CMS frame. No kine- 
matic cuts were applied to data. Green squares are A;, m 's from the 
GENBOD events. Gold inverted triangles are the SHD coefficients of 
Equation [24] for k=2. Black circles, blue stars and red triangles are 
SHD coefficients of the first, second and third terms, respectively, of 
the right side of Equation [25] Open circles are SHD coefficients of 
the right side of Equationl25l 



FIG. 10: (Color online) SHD coefficients for GENBOD -generated 
events consisting of 18 pions having average kinetic energy per par- 
ticle K = 0.9 GeV, as measured in the pair CMS frame. Only parti- 
cles with |r)| < 0.5 used in the correlation function. Green squares 
are A; m 's from the GENBOD events. Gold inverted triangles are the 
SHD coefficients of Equation [24] for k=2. Black circles, blue stars 
and red triangles are SHD coefficients of the first, second and third 
terms, respectively, of the right side of Equation[25] Open circles are 
SHD coefficients of the right side of Equationl25l 



We start with the case of interest- D = 3 spatial 
dimensions- and conserve 3 -momentum p. We implement 
energy conservation and on-shell constraints a bit later. 

We define JH 



f(Pi) 



d 3 N 



(6) 



as the single-particle momentum distribution unaffected by 
EMCICs. This may be considered the unmeasured "parent" 
distribution. Then, the k— particle distribution (k less than the 
total multiplicity AO including EMCICs is 



We implement total energy conservation = -y/s, by re- 
placing 8 3 (L/LiPi) — * § 4 (Y.'iLiPi P) m Equation[7] Here, 
P = ( t/s, ) is the total energy-momentum of the event, and 



P0,i 



pf + mj is the energy of the on-shell particle. 



fc{p\,-,Pk) = \\fiPi 



i(nU d3 Pjf(pj)) 53 $liPi) 



(7) 



Note the difference between numerator and denominator in 
the starting value of the index j on the product. 



We denote Lorentz-invariant distributions as 



f( Pl ) = 2E~=2E i f(p i ) 
dpf 



(8) 
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and rewrite Equation|7]as 



where 



fc(pi,...,Pk)=[Uf(Pi) 



S{Yf i . k+ ^~f{p J ))?> A ^Pi-P 




/( 




)f(p J ))s 4 (Ll l P.-P) 






[p)- m2 i) 


fiPj)) 


^ 4 (LI,p,-p) 



UfiPi 



I 


[u% k+i d 4 P jg (Pj) 


)&(Ll iP i-P) 


f(nU d *pj8(pj)) 


8 4 {LliPi-P) 



(9) 



fd 4 pg(p)-p n , 
fd 4 pg(p) 



i \ _ i ' d A pg{p)p hl pv 



(12) 



Finally, in Equation[lO] |B| denotes the determinant of the ma- 
trix B. Up to a factor of M, B is the inverse of the covariance 
matrix of the distribution g(p): 



B 



1 

M' 



-.b 



//V • 



(13) 



(14) 



We can now apply the CLT by recognizing the integral in 
the numerator in Equation [9] as the distribution of N — k mo- 
menta Ylj=k+iPj = P~uj=\Pj so tnat f° r "l ar ge enough" 
N — k,we find 



Thus, we arrive at an integral over four independent variables, 
in which the integrand function g(p) is "highly peaked" and 
with strong correlations in the 4-d p— space. 

According to Equation |9j the £-body momentum distribu- 
tion, including EMCICs, is the £-body distribution not af- 
fected by EMCICs- i.e. just an uncorrected product of 
single-particle distributions- multiplied by a "correction fac- 
tor" which enforces the EMCIC. The numerator of this factor 
counts the number of configurations in which the remaining 
N — k on-shell particles conspire to conserve total energy and 
momentum, and the denominator normalizes the distribution. 



B. Application of the Central Limit Theorem 

To arrive at a useful result, we argue along lines similar to 
those of 11351 l36l 13711 . The distribution of a large number M 
of uncorrected momenta W = Y4L1 Pi i s > by me Central Limit 
Theorem, a multivariate normal distribution 



M \ I M \ 

F M (W) = I [ l\d 4 P , g ( Pi ) J 8 4 I £ Pi -W (10) 




exp(--(W-(P"))^ v (W v -(P v )) 



Here, the average of the sum of 4-momenta is simply related 
to the single-particle average of the 4-momenta as 

M 

(P") = I(^)=M(^}, (11) 



fc(pu..,Pk)=\Uf{pi)J- 



N 



Fn(P) 




N-k 



(15) 



exp 



It is appropriate at this point to repeat the two approxima- 
tions we have employed up to now. The first assumption, al- 
ways important in using the CLT, is that N — k is sufficiently 
large; recall that N is the total multiplicity and k is the order of 
the correlation being calculated (k = 2 for two-particle corre- 
lations). Secondly, we have implicitly assumed that all parti- 
cles in the system are governed by the same single-particle 
distribution g(p). Strictly speaking, then, the system must 
consist of particles all of the same mass, and if there are sev- 
eral species with the same mass (say, K~ and K + ), they must 
furthermore have the same momentum distribution. This is at 
best an approximation for hadron or ion collisions, in which 
other particles contribute to the pion-dominated final state. 



C. Observable EMCIC effects 

Even the single-particle momentum distribution is affected 
by EMCICs: 



fciPi) = f(pt)- 



N 



N-l 



exp 



(16) 
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The product of such a single particle distribution forms the 
denominator of the ^-particle correlation function 



C{p\,...,Pk) 



fc(Pl,-,Pk) 



' _N_\ 2 

\ N—k ) 
,2k 



exp 



fc{pi)--fciPk) (^) 

w^n^-^))b^{p)-{p v ) 



(17) 



-l 



exp 



In this paper we concentrate on correlation functions in 
lout Aside and qi otl g, as is done in femtoscopic studies. How- 
ever, the two-particle correlation function in relative az- 
imuthal angle, which probes elliptic flow, may also contain 
EMCIC contributions through Equation[l7] These effects turn 
out to be small and are discussed in AppendixlAl 

To first order in 1 /N, the two-particle correlation function 
becomes 



C{pi,P2) = 1 



1 t " 



{p u ))b tN {pl-{P V ))- (18) 



The multivariate CLT used in Section lVBl accounts for cor- 
relations between vector components via the covariance ma- 
trix b~ (Eg. IT4b which has, in general, 10 nonvanishing ele- 
ments. The average vector P (Eq.fTTIi has in general 4 nonva- 
nishing elements. We now reduce these numbers significantly 
by considering the specific case of our interest. 

Firstly, we choose to work in the global center-of- 
momentum frame, so that 



</»"> = 5„,o<£) 



(19) 



As for the correlations, we are interested in signals gen- 
erated by EMCICs alone, not, for example, dynamical cor- 
relations due to flow. Neglecting elli ptic flow (azimuthal 
anisotropics in the parent distribution |5, 39]) implies 



(20) 



(b- l ) l2 = ( PxPy )=o. 



The same approach was adopted in earlier work H35LI361I38H . 
Similarly, we assume no dynamical correlations due to di- 
rected flow J39D, implying 



1.3 



'2.3 



= 0. 



(21) 



The on-shell constraint generates an unavoidable depen- 
dence between energy and 3 -momentum components. How- 
ever, in the CLT limit, only the second moment (covariance) 
comes into play, and this vanishes. For i ^ 0, 



(b-') . = (Ep i )-(E){ Pi ) = (Ep i ) 



(22) 



JdEjd 3 p-Eg(p)- Pi 
fdEfdip-g(p) 



0. 



In the last step, we recognize that p\ is an odd function of 
momentum, whereas E and g are even. 



In this scenario of interest, then, b is diagonal, and Equa- 
tions [T6lbecomes 



fc(Pi)=f(Pi)- 

1 

exp 



N 



, N ~ 1 
[Pi 



(23) 



Ply 



Pi, 



(E,-(E)) 2 



2(N-l)\(pl) (fi) (pi) (E^)-(Ey 



Similarly, EquationfTTlbecomes 

n , X fc{Ph-;Pk) 

C{Ph---,Pk) = -j 



exp 



-l 



fc{p\)---fc{Pk) 



1 ^=1 \ 



( —) 2 

(JL.) 2k 

(E 2 )-{E) 2 



(24) 



exp 



and Equation[T8lbecomes 



C(pi,pz) = 1- 

_j_ (r, Pl.T ■ P2.T 



(25) 



PU-P2,z , (E 1 -(E))(E 2 -(E)) 



(P 2 z) 



(£2) 



where we have taken (p 2 ) = (p 2 ) = (pj) /2 in the 
azimuthally-symmetric case of interest. In what follows and 
in Figures [BTTOl we shall refer to the first, second, and third 
terms within the parentheses of Equation |251 as the "pr" "pz" 
and "E" components, respectively. 

If we somehow know N, (pj), {pi), (E 2 ) and (E), we can 
calculate EMCICs using Equation[24] (See, however, the dis- 
cussion at the start of the next Section.) Better yet, if N is large 
enough, then we can use Equation [251 This is what is done in 
Figures |5|[T0l The open circles and orange inverted triangles 
represent the results of Equation [24] and Equation|25] respec- 
tively. The black circles, blue stars, and red triangles show 
the individual components of Equation l25l this decomposi- 
tion will be relevant when we discuss the "experimentalist's 
formula" in the next Section. 

Figures IBTTOl make clear that each of the three terms in 
Equation |251 produces non-trivial behavior of the A/ m 's. Also 
clear is the importance of not neglecting the energy term. We 
find also that the p z term affects Ai£, this may be surprising 
since A^^i quantifies the behavior of the correlation function 
in the "out-side" plane, while z is the "long" direction in the 
Bertsch-Pratt system. Clearly, EMCICs projected onto a 2- 
particle space are non-trivial objects. 

The first-order expansion (Equation (25) agrees well with 
the full expression (Equationl24b well for N > 10. Such mul- 
tiplicities are relevant for p + p measurements at RHIC (re- 
calling that N includes all particles, even unmeasured ones). 
We see also that the analytic calculations (open circles and in- 
verted triangles) approximate the results of the GENBOD sim- 
ulation (green squares), especially as the multiplicity and to- 
tal energy of the event increases; increasing agreement for 
large N and E tot is expected, given the approximations lead- 
ing to our analytic expressions. We observe also that the 
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analytically-calculated expressions respond identically to the 
kinematic cuts as does the simulation (c.f. Figures[9]and ITOb . 

Finally, the analytic calculations never reproduce exactly 
the simulations; we discuss this further in the next Section. 



VI. AN EXPERIMENTALIST'S FORMULA 

Even for large N and energy, the calculations do not exactly 
reproduce the EMCIC effects in the simulation. One reason 
for this may be found, in fact, in the definition of the aver- 
age values (e.g. (pf)) themselves. In Equation fT2l average 
quantities are calculated using the distribution f(p), which 
is not affected by EMCICs. Naturally, the only measurable 
distribution available to the experimentalist (even when GEN- 
BOD simulations serve as the "experiment") is f c (p)- 

Thus, it appears the experimentalist cannot plug her data 
into the equations[T2]and[25]to fully calculate EMCICs. How- 
ever, such an ambition would have been hopeless anyhow. Af- 
ter all, even the total multiplicity N (again, including photons 
etc) is rarely fully measured, and in principle N is a number 
of "primary" particles, a murky concept in itself. 

To the practicing femtoscopist, there is a natural solution. 
Having at hand (1) educated guesses for the quantities N, (E 2 ) 
etc, and (2) a physically-motivated functional form which 
connects these quantities to the correlations, one may perform 
a fit. Let us rewrite Equation [25] as 



C(pi,p 2 ) = 1 -Mi ■ {pi tT ■ P2,t} -M 2 - {pi, z - P2, z } (26) 

Ml 



where 



-M 3 • {Ei ■ E 2 } +M 4 ■ {Ei + E 2 } - 



2 1 



M 3 = 



N(p 2 T ) 
1 



N(pl) 



N((E 2 )-(E) 2 ) 



M A 



N((E 2 )-(E) 2 ) 



2V ( 27 > 



The notation {X} in Equationl26"l highlights the fact that X 
is a two-particle quantity which depends on pi and p 2 (or q, 
etc). From a practical point of view, X will be averaged over 
the same q bins as used for the correlation function. For in- 
finitesimally narrow g-bins, {X} = X. The binned functions 
{X} then automatically reflect the same event and particle se- 
lection as the correlation function. This involves nothing more 
than adding four more histograms to the several already being 
constructed by the experimentalist as she processes pairs in 
the data. 

Here, we should emphasize that, in Equation [26] p\, E\, 
P2, and E 2 should be calculated in the collision center-of- 
momentum (CCM) frame. The reason is that Equation [9] 
(hence Eqs. IT0ll27t assumes some fixed total energy and mo- 
mentum to be conserved. The event's total energy and mo- 
mentum (hence (E), (p), etc appearing in Eqs. rToETl i are 
fixed quantities in any given frame. In a pa/r-dependent frame 
(e.g. LCMS), the total energy and momentum of the event will 
fluctuate, pair-by-pair. Thus, while the correlation function 



may be binned in whatever frame one chooses, the momenta 
pi on the right side of Eqs. [911771 must be calculated in a pair- 
independent frame. In fact, starting with Equation[l9] we have 
chosen the CCM, for simplicity. 

The parameters M, defined in Equation [27] on the other 
hand, are global and independent of pi and p 2 . It is these 
which we will use as fit parameters. The task is then fast and 
straightforward; the EMCIC part of the correlation function 
C(q) is simply a weighted sum of four functions. Indeed, one 
may calculate coefficients as in Equation [4] for the four new 
functions. For example 



A l Z m(Q)= L {pi,,- jP 2,J(e,COS0 i ,^)-i'/,m(cOSe ( ^ ( ), 
bins ( 

(28) 

etc. Then, thanks to the linearity of Equation [26] and the 
orthonormality of F/ m 's, the measured A/ ,„'s themselves are 
similarly just weighted sums of harmonics 

A,,,„{Q) = 8/.0 • (1 -Ml/Mi) -Mi -Af T m (Q) (29) 



-M 2 -A<;% (Q) -M 3 -A?/' (Q) +M 4 -A): n r> (Q) . 

Treating Equation[29]as a fit, we have a few (say six, for / < 4) 
one-dimensional functions to fit with four adjustable weights. 
The number of degrees of freedom in this four-parameter fit 
remains high: ~ 300, for six A/,„,'s, each with 50 bins in \Q\. 

An example is shown in Figure [TT] where GENBOD - 
calculated correlation functions are fitted with the form of 
Equation[29] Not surprisingly, the minimization procedure re- 
turned fit parameters M, very close to the values calculated via 
Equation [27] However, exact agreement between the "best" 
parameter values returned by the fit, and those from Equa- 
tion [27] is not expected. This is because the large-Af approx- 
imation is only approximately valid and due to the fact that 
f{p) 7^ fc(p)> as discussed previously. Treating the Mi as ad- 
justable parameters leads to a slightly different weighting of 
the terms, and a slightly better fit to the data. 

Our original goal was not to understand EMCICs perse, but 
to extract the femtoscopic information from measured two- 
particle correlations. Assuming that the only non-femtoscopic 
correlations are EMCICs, one may simply add the femto- 
scopic terms OfemtoCPbf^) to the fitting function in Equa- 
tion [26] or [29] 



E+E) 



(30) 



1 -Mi ■ {pij-p 2 , T } -M 2 ■ {pi, z -p2, z } 



- M 3 -{Ei-E 2 }+M 4 -{Ei+E 2 }-^y 

Common femtoscopic fitting functions (e.g. Gaussian in the 
out-side-long space) usually contain ~ 5 parameters (e.g. N, 
X, Rj, i = o,s,l). In the imaging technique i40ll . one assumes 
the separation distribution is described by a sum of splines, 
rather than a Gaussian; here, too, there are usually 4-5 fit pa- 
rameters (spline weights). So, by including EMCIC effects, 
we have roughly doubled the number of fit parameters, rela- 
tive to a "traditional" fit which ignores them. This is a non- 
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FIG. 1 1 : (Color online) Green inverted triangles show the ir~ — 7I _ 
correlation function, in the pair rest frame, from 18-pion GENBOD - 
generated events. The black curve is a result of a fit with the "ex- 
perimentalist's formula" from Equation [29] Other curves represent 
the three component terms of the fit: Mi -Af^^Q) in the brown dot- 
ted line; M 2 -Af z m (Q) in the blue dashed line; M 3 -AfJ? (Q) +M A ■ 
Af^ t E (Q) in the red dash-dot line. See text for details. 



trivial increase in analysis complexity. However, we keep in 
mind two points. 

Firstly, the increased effort is simply necessary. EMCICs 
(and possibly other important non-femtoscopic correlations) 
are present and increasingly relevant at low multiplicity. One 
option is to ignore them, as has sometimes been done in early 
high-energy experiments. However, with the new high-quality 
data and desire for detailed understanding at RHIC, ignoring 
obvious features such as those presented in fiHl is clearly un- 
acceptable. Perhaps a slightly better option is to invent an 
ad-hoc functional form 12511 without a strong real physical ba- 
sis. We hope that the results here present a relatively painless 
and more reasonable, third option. 

Secondly, while the non-femtoscopic EMCICs are not con- 
fined to the large-Q region (an important point!), the femto- 
scopic correlations are confined to the small-Q region. There- 
fore, one hopes that the addition of four new parameters 
to the fit of the correlation function will not render the fit 
overly unwieldy. While we cannot expect complete block- 
diagonalization of the fit covariance matrix, one hopes that the 
M, are determined well enough at high Q that the femtoscopic 



fit parameters can be extracted at low Q. 



VII. NON-IDENTICAL PARTICLE CORRELATIONS 

For at least two reasons, it is important to turn attention to 
correlations between non-identical particles. 

First, it is natural to ask whether one can use other parti- 
cle combinations to "correct" for effects of EMCICs in, say, 
identical-pion correlation functions. After all, EMCICs are 
induced by global constraints on the entire event, not a spe- 
cific particle species. For example, various experiments have 
explored using (k + , K~ ) correlations to account for EMCICs 
in (z + ,n + ) correlation functions H! |H @] . 

Second, it is also important to know whether EMCICs 
could cloud the interpretation of correlations between non- 
identical particles. It is increasingly common to study asym- 
metries in the correlation functions of, say n — K pairs 14511 . 
interpreting such as a "shift" in the average point of emis- 
sion between the two particles l46ll . In the spherical harmonic 
decomposition, such shifts appear in the / = 1 moments (c.f. 
Appendix|B). We will find that EMCICs can indeed generate 
an asymmetry which might naively be considered proof of a 
femtoscopic shift. 

Here we discuss two effects- one immediately obvious and 
one more subtle- which are relevant for the above issues. The 
discussion is broken into three parts. Neglecting EMCICs and 
any other source of correlation at first, we briefly show the 
effects of two common resonances on correlations between 
oppositely-charged pions in a toy model. Thus calibrated, we 
use the more realistic and complex PYTHIA model to illustrate 
a non-trivial interplay between EMCICs and the resonances, 
which can mock up a femtoscopic asymmetry signal. Finally, 
we return to a toy model- now with non-identical particles and 
EMCICs, but without resonances or the several other sources 
of correlation present in PYTHIA - to make clear the mecha- 
nism behind the special effects EMCICs have on non-identical 
particle correlations. 



A. Effect of resonances 

First we consider the effect of resonances. To focus on ef- 
fects other than global EMCICs we use a toy model in which 
only ten identical resonances per event are generated and no 
other particles. The momentum of each resonance is gener- 
ated from a thermal distribution; energy and momentum are 
conserved for each decay separately, but not globally for the 
whole event. 

Figure [12] shows the spherical harmonic moments of 
(7t + ,7t~) correlation functions for events including co (blue 
squares) and p resonances (red triangles). As seen, even with- 
out considering EMCICs, the correlations among particles 
coming from resonance decays produce non-trivial structure. 
In this case, one cannot simply divide (tt + , jt + ) correlation 
function by (k + ,k~) to "remove" EMCICs. 
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FIG. 12: (Color online) (7t + ,Jt~) correlation functions calculated in 
the LCMS frame for events including co (blue squares) and p (red 
triangles) resonances . 
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FIG. 13: (Color online) SHD moments of (lt + , Jt~) correlation func- 
tion from p + p and p + p collisions at 200 GeV calculated from 
PYTHIA events. 



B. Entrance channel asymmetries 

In addition to the correlation between daughters of reso- 
nance decays (c.f. Fig. [I2l . there is a more subtle effect to 
consider. This happens when the two particles have different 
inclusive momentum distributions and energy and momentum 
are globally conserved. Under these conditions non-identical 
particle correlations exhibit structure absent in identical parti- 
cle correlations. 

FigurefHlshows PYTHIA J47t] calculations of (jl + ,7r.~) cor- 
relations for p + p and p + p collisions at 200 GeV. In addi- 
tion to obvious correlations between daughters of resonance 
decays (K®, CO, p), we see additional structure. We focus on 
the structure in the / = 1 moments. In general, such moments 
need not vanish for correlations of non-identical particles, as 
discussed in AppendixlBl 

Correlations between sibling daughters of p and co reso- 
nance decay do not generate I = 1 moments, as seen in Sec- 
tion [VTlA] However, pions which are daughters of these de- 
cays will in general have a different single-particle momentum 
distribution than pions from other sources in the event. If the 
fraction of pions from resonance decay, as a function of pion 
momentum, is different for K + and 7t~, then the single-particle 
distributions of positive and negative pions will be different. 
We argue below that it is this difference in single-particle dis- 



tributions which is the key to the non-vanishing / = 1 mo- 
ments; that this difference may arise from resonances in the 
case at hand is irrelevant. 

In the p + p collisions, the fraction of n + coming from any 
given source (e.g. p-decay) must be identical to that of n~, for 
a given value of pj. Thus the pj distribution of K~ must be 
identical to that of K + . However, the rapidity distributions will 
be mirror images of each other. Thus, any asymmetry in K~ — 
K + correlations from p + p collisions will be associated with 
<7iong! and will appear in A^o, as seen in Figure [13] Similarly, 
the vanishing (non-vanishing) moment A\q (Ai i) for p + p 
collisions reflects the fact the single-particle distributions will 
show no asymmetry in rapidity, but may differ as a function 
of pj. 

Since single-particle distributions divide out of a correla- 
tion function, a difference between K + and K~ momentum dis- 
tributions, by itself, cannot generate a signal in A/, m 's. Rather, 
a global correlation, coupled with this difference, generates 
the signal. We discuss this further below. 



C. A simpler case 

In Section IVII B| we argued that the small difference 
in single-particle momentum distributions between positive 
and negative pions produced by PYTHIA , coupled with 
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FIG. 14: (Color online) Green squares show the % — p correlation 
function, in the pair rest frame, from GENBOD -generated events. 
Gold inverted triangles are the SHD coefficients of Equation 1241 for 
k=2. Black circles, blue stars and red triangles are SHD coefficients 
of the first, second and third terms, respectively, of the right side of 
Equation [25] Open circles are SHD coefficients of the right side of 
Equationl25l 



global conservation laws, generated non-trivial EMCICs in 
the non-identical particle correlations. However, PYTHIA con- 
tains many non-EMCIC sources of correlations, related to 
string fragmentation and other processes, which might be 
flavor/isospin-dependent. To make clearer our argument, we 
here show a simple GENBOD simulation, containing both pi- 
ons and protons, but no explicit correlations between them 
such as a A resonance. Due at least to their different masses, 
/proton 7^ /pion is guaranteed. 

Figures [14] shows the K — p correlation function. Since 
the underlying single-particle proton and pion distributions 
are isotropic, Aio (sensitive to shape elongation in C(q) in 
g/relative to transverse components) is expected to vanish. 
Ai i is finite, however, due to differences in pj distributions. 
Since there is no other source of correlation in the simulation, 
this obviously is an EMCIC. 

From Figure [14] it is also clear that neither Equation [24] 
nor its first-order expansion[25]fully describes the correlation 
function. This is due to the fact that our formalism is built on 
the assumption that all particles in the system follow the same 
parent distribution, as pointed out after Equation[T2l 



VIII. SUMMARY 

To truly claim an understanding of the bulk nature of matter 
at RHIC and the LHC, a detailed picture of the dynamically- 
generated geometric substructure of the system created in 
heavy ion collisions is needed. It is believed that this substruc- 
ture, and the matter itself, is dominated by strong collective 
flow. The most direct measure of this flow is a measurement 
of the space-momentum correlation (e.g. Rimj)) it gener- 
ates. The physics of this large system, and the signals it gen- 
erates, should be compared to the physics dominating p + p 
collisions, as is increasingly important in high-pr studies at 
RHIC. For small systems, however, non-femtoscopic effects 
contribute significantly to the correlation function, clouding 
the extraction and interpretation of femtoscopic ones. 

We have discussed a spherical harmonic representation of 
the correlation function which clearly separates components 
of the three-dimensional shape measured in modern experi- 
ments. This representation is maximally efficient, inasmuch 
as only a few one-dimensional plots need be examined to ex- 
tract full 3-D shape information. The relevant number of such 
plots is further reduced due to symmetry conditions, discussed 
in detail in AppendixlB] 

EMCICs, correlations generated by kinematic conservation 
laws, are surely present and increasingly relevant as the event 
multiplicity is reduced. Using the code GENBOD to study cor- 
relation functions solely driven by EMCICs, we found highly 
non-trivial three-dimensional structures strongly influenced 
by event characteristics (multiplicity and energy) and kine- 
matic particle selection. 

We extended the work of Danielewicz, Ollitrault and 
Borghini to include four- momentum conservation and applied 
it to correlation functions commonly used in femtoscopy. We 
found structures associated individually with the conservation 
of the four-momentum components, which interfere in non- 
trivial ways. Comparison of the analytic EMCIC calculations 
with the GENBOD simulation gave confidence that the approx- 
imations (e.g. "large" multiplicity AO entering into the calcu- 
lation were sufficiently valid, at least for multiplicities consid- 
ered here. We further showed that the full EMCIC calculation 
can safely be replaced with a first-order expansion in 1/N. 

Based on this first-order expansion, we developed a practi- 
cal, straight-forward "experimentalist's formula" to generate 
histograms from the data which are later used in a generalized 
fit to the measured correlation function, including EMCICs 
and femtoscopic correlations. The degree to which this func- 
tional form fully describes measured experimental correlation 
functions has not been discussed and will need exploration on 
a case-by-case basis. 

There is strong interest in correlations between non- 
identical particles, for two reasons. Firstly, sometimes K~ — 
K~ correlations are divided by n + — K~ correlations in an at- 
tempt to "divide out" EMCICs. (In such a procedure, reso- 
nance regions are avoided, naturally.) We discussed potential 
problems with such an approach, related to entrance-channel 
asymmetries coupled with EMCICs. Secondly, 3-D asymme- 
tries (sometimes quantified as "double ratios", e.g. l45ll ) in the 
correlation function for different-mass particles (e.g. K — K) 
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FIG. 15: (Color online) V2(pr) for different event multiplicities. See 
text for details. 



are often interpreted in terms of dynamically-generated differ- 
ences in the average space-time emission point between the 
two particles. Using a very simple example, we discussed that 
EMCICs might significantly cloud such an interpretation. 

The huge systematics of results and interest in femtoscopy 
in heavy ion collisions is renewing similar interest in the 
space-time signals from p + p collisions. Direct comparisons 
between the two systems are now possible at RHIC and have 
already produced intriguing preliminary results. Very soon, 
p + p collisions will be measured in the LHC experiments, 
and the heavy ion experimentalists will be eager to apply their 
tools. The femtoscopic tool is one of the best in the box- so 
long as we keep it sufficiently calibrated with respect to non- 
femtoscopic effects increasingly relevant in small systems. 
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APPENDIX A: EMCIC EFFECTS ON v 2 

Since EMCICs can produce a structure in the correlation 
function even in the absence of femtoscopic correlations, it is 
worthwhile to check analytically and then confirm with sim- 
ulations whether vi- a common measure of collective elliptic 
flow H]- may be affected by EMCICs. 



cos (mA^j cos (nA(j)) c/A(j) = 5 mj „Jt, 



(Al) 



where for v%, n = 2, 

This means that in the absence of flow all EMCIC terms 
vanish except for the ones that exhibit cos (2A(j)) dependence 
of A<|). For example in the first order expansion of EMCICs 
(see Eq. l24t there is a term pt,\Pt,\ ~ cos (A(|)). This term 
gives no contribution to V2, nor do any other terms from 1 /N 
expansion. The first term that gives a non-zero contribution to 
V2 (means, goes like cos (2A(j)) ) is the second order expansion 

term in pj that is proportional to (Pt,iPt,i) 2 ~ c °s 2 (A<|>) ~ 
cos (2A(|)). This term (as well as a few other terms in higher- 
order 1 /N expansion) will give a non-zero contribution to V2. 
In our GENBOD simulations we do not have a flow so we can 
study the magnitude of the EMCIC effects on v 2 measure- 
ments. Such results are presented on Figure Q3] where we plot 
vi vs pj for three different event multiplicities while the free 
kinetic energy per particle is fixed (K = 0.9 GeV). 

As seen, the magnitude of a non-flow contribution to V2 
from EMCICs is getting smaller with increasing multiplicity 
and even for low-multiplicity events the magnitude is of order 
of a few per-mile for large pj. From this dependence we can 
predict that this effect will be so small in heavy ion collisions 
that it can be simply neglected. 



APPENDIX B: SYMMETRY CONSIDERATIONS 

The spherical harmonic decomposition representation, in 
which three-dimensional correlation functions are represented 
by several one-dimensional moments, A/ m , efficiently con- 
denses the shape information. A much greater increase in ef- 
ficiency comes, however, with the realization that many A; m 's 
must vanish by symmetry, depending on the cuts and condi- 
tions of the analysis. Besides reducing information by signif- 
icant factors, this realization also provides diagnostic power- 
non-physical artifacts often appear in A/,,„'s which do not van- 
ish when they should. Digging out such effects in the tradi- 
tional three-dimensional Cartesian representation can be quite 
difficult. 

In the most general case, the 3-D correlation function may 
have any shape, with no symmetry constraints. In this case, 
none of the A/ m 's need vanish. Usually, however, an analysis 
is less than fully general, and symmetry consequences then 
arise. 

In particular, we will consider four common conditions 
used in practice: 

[A] , One measures correlations between identical particles 

[B] . The measurement covers a symmetric rapidity region 

about y=0 and the collision is between identical ions 
(e.g. Au+Au rather than Au+Si) 

[C] . The measurement is integrated over reaction-plane an- 

gle 
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# 


transformation 


Yl „,consequence 





(9o,9s,9l) -> (+9o,+9s,+9l) 




1 


(9o,9s,9l) -> (+9o,+q s ,-9l) 


!/,,„ ^(-l) /+ '"^, m 


2 


(9o,9s,9l) -> (+«o>-?s,+?/) 




3 


(q ,q s ,qi) -* (+q ,-q s ,-qi) 


¥,, m ^(-l) l+m Y,, m * 


4 


(9o,9s,9l) -> (-9o, 


Y l,m — * ( — 1 


5 


i9o,9s,9l) -* {-9o,+9s,-9l) 




6 


(9o,9s,9l) -* {-9o,-9s,+9l) 


*z,m^ (-l) ffl i7. m 


7 


(9o,9s,9l) -* {~9o,-9s,-9l) 


*Z,m^(-l)%» 



TABLE I: The possible transformations (numbered in the left col- 
umn) in which the signs of q components flip, and the effect of the 
transformation on the Y\ m 's. 

[D]. The measurement might be correlated with the second- 
order reaction-plane, but the first-order reaction-plane is 
not known. In other words, the direction of the impact 
parameter is known at best only modulo 71. 

Our strategy begins by identifying transformations in rela- 
tive momentum q under which the measured correlation must 
be invariant. As an example, since the overall sign of q is 
meaningless when discussing pairs of identical particles (con- 
dition |[A]) , C(q ,q ,qi) = C(—q ,—q o ,—qi), or, in spherical 
coordinates, C(2,cos8,(j)) = C(Q, — cos8,(|) — Jt). 

We then use a symmetry of the spherical harmonics, here 
Yi Jn (cos9,(j)) = (-l) f y i)m (-cos6,<|)+7C) to find 

A Lm {Q) = J'^dtyJ 1 JcoseC(e,cos9,4))y/, m (cos9,<|>) 
1 r 2n r l 

= — / d§ / dcosQC(Q, — cos9,(|) — 7t)7; m (cos8,(|)) 
4% Jo J -i 

i r n /•- 1 

= — / d& / (-dcose)C(<2,cose,(b)F/ m (-cose ! d)+Jt) 
47W-7I Ji 

= — / dty rfcos9C(e,cos9,<|))(-l)'y, m (cose,(|))) 

471 Jo i-1 

= (-l)%, m (<2). (Bl) 

Thus, all odd-/ moments A/ „,must vanish, for correlations be- 
tween identical particles. 

The same type of reasoning is used below, in identifying 
symmetry constraints for various combinations of analysis 
conditions. 



1. q transformations and Y t m response 

Table J] lists all combinations in which one or more of the 
components of q can change sign. For later reference, the 
transformations are numbered 0. . . 7, according to a binary 
scheme. The effect of the transformation on the spherical har- 
monics appears in the last column of the Table. 

Transformation (0), of course, is the trivial identity trans- 
formation, under which any correlation function is invariant, 
and which imposes no symmetry constraint. We include it in 
the Table only for completeness, and do not discuss it further. 
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TABLE II: Symmetry consequences of analysis conditions. The 
left four columns show various combinations of analysis cuts and 
conditions, identified [A]-[D] as discussed in the beginning of this 
Appendix. (Note that condition | [C] | implies condition |[D]| this is 
indicated by the symbol (^/) in column f[D]| ) The middle seven 
columns indicate the consequent invariance symmetries of the cor- 
relation function according to the numbering scheme of TableU The 
right-most column indicates which, if any, spherical harmonic mo- 
ments of the correlation function must vanish. 

2. Restrictions, invariants, and consequences on A/ m 's 

Under which of the transformations in Table|I]does the cor- 
relation function remain invariant? Since identical-particle 
correlations are more common than correlations between non- 
identical particles, there will be a greater familiarity with the 
symmetries of the former. Thus, we begin with this more 
familiar case and then discuss non-identical particle correla- 
tions. 



a. Correlations between identical particles 

To systematically identify those transformations in Table U 
which leave a correlation function invariant, it helps to have 
a concrete functional form to discuss. For identical pions, the 
correlation function is often parameterized as a Gaussian with 
six "radius" parameters, 

C(q a , q , qi ) = 1 + A, • exp {-R 2 q 2 - R]q] - Rjqj - (B2) 
-2R 2 0S q q s - 2R 2 ol q q t - 2R 2 sl q s qi) . 

While measured correlation functions often have non- 
Gaussian features not captured by this parameterization, the 
form given in Equation lB2l contains the generic and most gen- 
eral symmetries of all correlation functions using identical 
particles. Thus, we use this familiar example to focus the dis- 
cussion. The six parameters of form IB2 1 describe an ellipsoid 
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described by three axis lengths, and rotated by three Euler an- 
gles in g-space. Measured examples are shown and discussed 
in HI. 

Clearly, the form of Equation lB2l is invariant under transfor- 
mation (7), as discussed earlier. Invariance under any trans- 
formations (l)-(6) requires that one or more of the "radius" 



parameters Rfj vanish. In general, none of them do II48L l49tl . 
even when considering a region symmetric about mi drap idity 
in a collision between identical ions (condition |[B]| ) 15911 . 

If the measurement is integrated over reaction plane angle, 
then the "side" direction has no relevant sign, R 2 S = = 0, 
and the correlation function is invariant under transformation 



(2). While R 2 o[ netd not vanish 15011 , the correlation function 
is unchanged if ^ and ^/change sign together (transformation 
(5)). 

Further constraining the measurement to a symmetric re- 
gion about midrapidity implies also that R 2 , vanish, and the 
correlation function is then invariant under all transformations 
(0)-(7). This is the most common set of measurement condi- 
tions. 

At high energies, it is common only to determine the 
second-order reaction plane. This corresponds to condi- 
tion |[D]| If the measurement is performed at midrapidity 
(condition |[B]| i, then Rg S is the only non-vanishing cross-term 
radius, so the correlation function is invariant under transfor- 
mation (1). Away from midrapidity, R^need not vanish, so 
(7) is again the only remaining transformation leaving C(q) 
invariant. 



b. Correlations between non-identical particles 

Correlations between non-identical particles are no longer 
invariant under transformation (7), as they may depend on 
odd-power terms of the components of q. In the case of fem- 
toscopic correlations, the strengths of these odd powers probe 
asymmetries in the average emission point between the two 
particle species 114611 . 

From a symmetry standpoint, the correlation function will 
be characterized by nine parameters, rather than the six "HBT 
radii" of Equation IB2I In the simple case that C(q ,q ,qi) 
would be Gaussian, these new parameters might represent the 
offset from the origin of the ellipsoid in g-space. 

In the absence of any cuts- or if only the midrapidity condi- 
tion |[B]| is applied 10, a ll nme parameters may take any value, 
and there are no required invariances or symmetry constraints. 
If the reaction-plane is integrated over (condition |[C]| l then 
C(q () ,q ,qi) may remain sensitive to the sign of q (reflect- 
ing, for example, a different average time of emission between 
the particles 14611 ) and q\ (reflecting the difference in emission 
point in the beam direction, for analyses away from midrapid- 
ity), but not q s , since an angle-averaged physical source must 
be symmetric with respect to the beam axis. 

Unlike the case in which it is the sole condition, if the 
midrapidity condition |[B]| is imposed together with condi- 
tion |[C]| then it does have an effect. In particular, a depen- 
dence on the sign of ^/vanishes. 



If condition |[C] | is relaxed to condition |[D]| (i.e. the anal- 
ysis is sensitive to the second-order reaction plane), then the 
sign of q s may matter. This is because the sign of q always 
affects correlations between non-identical particles and, as in 
identical particle correlations in which R 2 S may be finite, so 
the sign of q q s may separately matter. Thus, imposition of 
|[Dl| alone implies no symmetry constraints. 

APPENDIX C: FINITE BINNING EFFECTS 

Equation|4]defines the harmonic moments in terms of a con- 
tinuous correlation function. Most experimentally-measured 
correlation functions are constructed via histograms with dis- 
crete, finite bins. For decomposition into spherical harmonics, 
a natural choice would be to use bins in Q, cos0 and (j) (c.f. 
Eq.[3]l. Here, we will find an approximate expression, analo- 
gous to Equation|4] for the harmonic moments in terms of the 
discretized correlation function. 

We denote the fixed bin sizes in the angular coordinates as 
A cos e and A^. Binning in Q is unimportant here, since Q is 
carried as an explicit argument in both C and A/ ,„. The binned 
correlation function (denoted with superscript A) is related to 
the continuous one as 



C A (e,cos9,^,) 

I / -cos9,+Acos0/2 



AaA, 



cos9 Jtyi— A^/2 Jcos9,— Acos9/2 



J(cos9)C(e,cos9,(|)) 



+/' 



l^—L E A lW (Q) 

^0^COS9 ;/_n _,/__// 



X 



l'=0m'=-l' 

/■cos9,+A cose /2 

d$ / £f(cose)y; m ,(cose,<|>) 

,-V2 icos9,-A cos9 /2 
_ == +/' 

471 E E A />'(2)- F />'( A <|)Ao S e,cose i )x 

l'=0m' = -l' 

y; m ,(cose,-,cK). (ci) 



Here, 

sin(mA<j,/2) 
Fi',m' j Acos e , cos 6,- J = — ^ — x 



^cos 

eP/' ; „,'(cos9,) . / cos9,-A cos9 /2 



cos9,+A cos9 /2 



(C2) 



rf(cos0)P//, m /(cos0) 



is the term which includes the finite binning effects. 

Assuming that A/ m 's vanish for l,m greater than the sam- 
pling Nyquist frequency, by the sampling theorem lf5ll l52ll . 
the A/ m 's are completely determined by C A . In fact, if F[ m 
were independent of cos0,, then we would have 



A, m (Q) 



4^008 9 



Pi,m (AosejAfr) \/4rc 
^C A (e,cos0,^,)y / ,„(cos0,^,) : 

bins / 



where the summation is over all bins of cos and (j), for a given 

Q. 
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However, F\ m does depend on cos 0,, so the above equation 
does not strictly hold. Nevertheless, we find, numerically, that 
an excellent approximation is 

A, (Q) ~ A( ^ Acos9 y C A (Q,cos&i,§i) Y Um (cos 9j , 0,- ) 

V47I Ft >m (Acqs 6,^,0086;) 

(C3) 

For any given measurement, one may check the validity of 



this approximation by plugging the result of Equation lC3l into 
the expression on the last line of Equation |C II To the extent 
that it returns the measured correlation function C A , the A/ m 's 
returned by Equation IC3I are correctly extracted. If there are 
deviations, the correct A/,,„'s can be found by iterative tech- 
niques. 

Other methods to remove binning effects have also been 
proposed |[53ll . 
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